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1. Introduction 

1.1 Motivation 

Twisted mass lattice QCD (tmLQCD) offers an efficient mechanism for eliminating unphys- 
ical zero modes[|l|] and removing 0(a) errors from Wilson simulations [g]. An advantage of this 
formulation is that it allows simulations to be done at relatively light quark masses. Good results 
can be obtained for pions[^, ^ but the correlation functions for other hadrons tend to have sig- 
nificant statistical fluctuations when the quarks are light. In particular our previous exploration of 
the octet and decuplet baryons indicated that simple local operators are inadequate at small quark 
masses [S]. Here, we study whether the application of sink smearing can be useful in improving 
baryon signals and reducing statistical errors. 

1.2 Action 

We use the Wilson gauge action with a tmLQCD fermion action^] 

S f [Y,xjf,U] =a A Y l W) (^ + ^75 + 7-V ± -^£v;V v ^) ¥ (x), (1.1) 

where 

y T = (u,d,c,s), (1.2) 

M = diag(Mi,Mi,M c ,M s ), (1-3) 
PL = diag(n,,-n h n c ,-n s ). (1.4) 

Forward (V), backward (V*) and symmetric (V^) derivatives are standard, and the c field is irrele- 
vant for our quenched simulations. 

1.3 Simulation details 

Simulation parameters are given in Table [j] These parameters correspond to maximal twist 
using the parity definition^, [7|, §]. In the following, multi-state fits are used to extract ground state 
masses. 



2. Smearing Method 

Gaussian quark smearing is used at the sink[^], 

q smr (x) = (l + aA) na q(x), (2.1) 

&q(x)= £ [Uj(x)q(x + j)-q(x)]. (2.2) 

y=±l,±2,±3 



with stout spatial link variables [10] 



1/(2) _>..._»1/(»p) ) (2.3) 
U ( j n+1) (x) = exp (ip® ( f(x)) U { j"\x). (2.4) 
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The four smearing parameters were optimized for the reduction of excited state contributions at 
(5=6.0, leading to a=0.15, n a =64, p=0.15 and n p =\6. The link smearing parameters were opti- 
mized first using the static quark anti-quark potential (see Figure [I] and compare with [11]). A 



comparison between the effect of quark smearing and link smearing in reducing excited state con- 
tributions at our lightest and our heaviest masses is shown in Figure |2| 

3. Results 

3.1 Vector Mesons 



In [|12|], results for the pseudoscalar and vector mesons with local operators were reported. A 
mass splitting between charged and neutral kaons was observed but, because of the large errors, no 
similar splitting between the neutral and charged vectors was found. 

For any quark flavor q, the operators qyjq and qOj^q both couple to the neutral vector meson, 
but the latter can mix with the axial vector away from maximal twist. In the present work, sink 
smearing was implemented for both of these operators, and the charged-neutral splitting was still 
insignificant at these lattice spacings. Representative effective mass plots are shown in the left 
panel of Figure ||. 

Continuum extrapolations were performed linearly in a 2 for each quark mass. Subsequently, 
a simple linear fit in quark mass was used to approximate a chiral extrapolation, as shown in the 
right panel of Figure ||. Masses for the p and K* that were obtained in this way are 0(10%) above 
experiment, as expected for a quenched simulation. 

3.2 Spin 1/2 baryons 

Source smearing allows for a more precise determination of ground state octet baryon masses. 
(For our previous exploration with local operators, see From the left panel of Figure ^, one 
sees the good scaling of the £~ mass with lattice spacing, while the Z + has a significant linear 



j3 


#sites 


#configs 


aM 


a\x 


physical mass 


5.85 


20 3 x 40 


594 


-0.8965 


0.0376 


~ m s 








-0.9071 


0.0188 


~ m s /2 








-0.9110 


0.01252 


~m s /3 








-0.9150 


0.00627 


~ m s /6 


6.0 


20 3 x 48 


600 


-0.8110 


0.030 


~ m s 








-0.8170 


0.015 


~ m s /2 








-0.8195 


0.010 


~ m s /3 








-0.8210 


0.005 


~ m s /6 


6.2 


28 3 x 56 


200 


-0.7337 


0.021649 


~ m s 








-0.7367 


0.010825 


~ m s /2 








-0.7378 


0.007216 


~ m s /3 








-0.7389 


0.003608 


~ m s /6 



Table 1: Simulation parameters. 
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Figure 1: Effective mass for j3 = 6.0 at temporal separation t=a corresponding to the static quark-antiquark 
potential at spatial separation r=5a, as a function of the stout link smearing parameters p and n p . The fitted 
ground state shown with no smearing used a two-exponential fit. 



dependence on dr. Since our tmLQCD corresponds to degenerate up and down quark masses, the 
observed splitting among E masses at j3 = 5.85 is a twist artifact. For j8 > 6.0, this artifact is not 
statistically significant. A similar twist artifact was observed among the E baryons. 

Like the vector mesons, the quenched baryon masses obtained from simple linear fits in the 
quark mass are systematically larger than experiment, as evidenced by the right panel of Figure |j. 
The ordering of N, A, £ and S masses is correct. 

3.3 Spin 3/2 baryons 

Even with smearing, quarks much lighter than m s make fitting difficult. In Figure |5[ the best 
case (£2 has only strange valence quarks) and the worst case (A has only light valence quarks) are 
shown. As expected, the quenched Q. is heavier than experiment. The A data are noisier than one 
would like; clearly sink smearing and/or additional operators should be explored in this context. 
Hints of twist artifacts can be seen at j8 = 5.85, but the over-all poor quality of A(1232) data 
underscores the importance of future improvements like group-theoretical operators and source 
smearing. 



4. Discussion 



Sink smearing was effective in reducing excited state contributions. However, this did not 
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Figure 2: Effective mass plot for the charged vector at j3 =6.0 with the (a) heaviest and (b) lightest quark 
masses. Unsmeared results (black), quark smeared results (blue), and quark and link smeared results (red) 
are shown. 



* * i 



. y (local) 
a o j4 (local) 



• y. (smeared sink) 
a Oj 4 (smeared sink) 



1 1 1 ' 


i 1 i 1 


i 1 i 1 


: >' 




" «i> 






- p 
i.i. 


i.i. 


i.i. 



0.2 0.4 0.6 0.8 1 

quark mass in units of m 



Figure 3: The left panel shows effective mass plots for the neutral vector meson with j3 = 5.85 and quark 
mass = m s /2 for two local operators. The right panel shows simple linear fits to vector meson masses after 
continuum extrapolation. 
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Figure 4: The left panel shows scaling of the E + and E~ masses for three light quark masses, as well as 
linear extrapolations in a 2 . Tiny horizontal offsets are for clarity. The right panel shows simple linear fits to 
the octet baryon masses, after continuum extrapolations. 




Figure 5: Scaling of the Q. and A masses. Linear fits to the A data are illustrative; note that j3 = 6.0 and 6.2 
scale within statistics. Tiny horizontal offsets are for clarity. 



always lead to a more precise prediction for the ground state mass, compared with multi-state fits 
to unsmeared data. Nevertheless, some improvement in reducing systematic errors arising due to 
multi-exponential fits may have been realized. 

Reasonably precise results for vector meson and spin 1/2 baryon masses were obtained, and 
indications of twisting artifacts were identified on our coarsest lattice. Even with smearing, the 
spin 3/2 data were of low quality for m u ^ < m s /2. Negative parity baryons and excited states were 
included in the multi-state fits, but were not useful for studies of the corresponding physical states 
and are not shown here. 

Sink smearing is useful but not sufficient to allow precise studies of the baryon spectrum with 
tmLQCD. In future studies, source smearing should be considered. As well, operators constructed 
according to the lattice symmetries are being developed for tmLQCD [ 13]. 
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